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ABSTRACT 


A  method  of  analysis  of  an  integrated  dielectric  slab  waveguide  and 
wedge  radiator  is  presented.  Plane  wave  constituents  in  the  slab  are 
assumed  to  excite  the  dielectric  wedge.  The  solution  in  the  wedge  is  pos¬ 
tulated  to  be  a  finite,  continuous  spectrum  of  plane  waves.  Reflections 
of  these  waves  from  the  wedge  boundaries  are  taken  into  account  by  intro¬ 
ducing  angularly  dependent  Fresnel  reflection  coefficients  into  the  inte¬ 
grands  of  the  plane  wave  integral  representations.  These  integrals  reduce 
asymptotically  to  the  dielectric  slab  fields  in  the  limit  of  small  wedge 
apex  angle  at  the  interface  between  the  slab  and  the  wedge. 

Using  numerical  techniques,  the  TE  surface  fields  and  the  field  ra¬ 
diation  patterns  are  determined  for  different  dielectric  materials,  wedge 
lengths  and  slab  widths.  It  is  observed  that  by  decreasing  either  the  rela¬ 
tive  permittivity  of  the  dielectric  or  the  slab  width  or  by  increasing  the 
wedge  length  a  more  directive  antenna  pattern  results. 
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INTRODUCTION 


1 . 

The  current  availability  of  low  loss  silicon  has  provided  the  capabil¬ 
ity  of  constructing  miniature  integrated  waveguide  and  antenna  devices.  Such 
devices  are  of  interest  to  U.S.  Army  scientists  for  use  in  the  millimeter 
wave  transmitters  and  receivers  being  developed  to  improve  data  transmission 
rates  [I],  Since  the  antenna  structures  under  consideration  are  constructed 
from  dielectric  materials  (silicon),  they  are  tapered  to  improve  their  radia¬ 
tion  characteristics  (increased  pattern  directivity  and  lower  side  lobe  levels) 
over  a  wide  frequency  band  [1-4).  Although  rigorous  theoretical  approaches  to 
study  tapered  dielectric  structures  in  principle  are  available  [12,13],  the 
solution  to  the  dielectric  wedge  problem  has  not  been  worked  out  in  detail. 

Tapered  dielectric  geometries  have  been  analyzed  using  local  mode  type 
theories  [6-11].  In  these  theories,  the  guided  wave  characteristics  of  a 
Capered  structure  are  determined  from  the  corresponding  waveguide  characteris¬ 
tics  of  the  untapered  structure.  A  first  order  solution  based  on  Shevchenko's 
rigorous  approach  [12]  yields  local  mode  results.  3oth  Marcuse  [8,9]  and,  to 
first  order,  Shevchenko  [12]  obtained  approximate  solutions  for  the  field  along 
a  tapered  slab  waveguide  joining  two  dielectric  slabs  of  different  heights  in 
terms  of  a  local  surface  wave  by  using  a  complete  set  of  transverse  mode  func¬ 
tions  and  presented  different  methods  for  determining  the  expansion  coefficient 
Balling  [6,7,14]  examined  the  near  and  far  fields  of  a  line  source  imbedded  in 
an  infinite  two-dimensional  dielectric  wedge.  He  compared  a  VKB  solution  to 
a  more  accurate  one  based  on  plane  wave  superposition  integrals  which  allowed 
for  ray  optical  and  lateral  ray  interpretations.  Methods  using  ray  optics  to 
determine  asymptotic  (ray)  modal  fields  in  non-radiating,  non-uniform  or 
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capered  geometries  iwith  boundary  conditions  which  allow  solutions  to  be  ob¬ 


tained  bv  separation  of  variable  techniques)  have  been  thoroughly  studied  by 
Maurer  and  Felsen  [15-13],  Their  approach,  however,  does  not  apply  directly 
to  the  dielectric  wedge.  Bates  [5]  confirmed  the  conclusion  of  Maurer  and 
Felsen  chat  source-free  modal  solutions  for  dielectric  wedge  structures  can¬ 
not  be  constructed.  The  rigorous.  Full  Wave  approach  of  Bahar  [13]  can  be 
used  to  study  the  wedge  problem,  but  it  is  extremely  difficult.  Shevchenko's 
rigorous  approach  [12],  which  is  similar  to  Bahar 1 s, would  require  higher  order 
terms  to  be  applicable  to  the  wedge  problem.  An  excellent  survey  of  work  done  on 
the  dielectric  wedge  problem  can  be  found  in  Balling's  dissertation  [6]. 

The  method  of  analysis  presented  here,  which  provides  meaningful  physi¬ 
cal  insight  into  the  wave  processes  taking  place  in  a  dielectric  slab-wedge 
integrated  structure,  is  based  on  the  work  of  Maurer  and  Siwiak  [19,20].  In 
that  theory,  a  finite  wedge  segment  is  excited  by  a  single  slab  mode.  The 
field  in  the  dielectric  wedge  region  is  inferred  from  an  integral  representa¬ 
tion  of  a  Hankel  function  whose  asymptotic  form  closely  matches  the  ray  struc¬ 
ture  of  the  slab  mode  in  the  common  region  between  the  slab  and  wedge.  Inter¬ 
preting  this  asymptotic  solution  in  terms  of  rays  leads  to  a  ray-tracing  proce¬ 
dure  whereby  rays  are  multiply  reflected  from  the  wedge  walls.  The  amplitude 
and  phase  of  the  reflected  ray  fields  are  assumed  to  be  altered  by  Fresnel 
reflection  coef f icients.  From  this  ray  field  an  approximate  expression  for 
the  actual  field  on  the  wedge  surface  is  then  formulated  in  terms  of  a  series 
of  integrals  whose  integrands  contain  saddle  points  or  stationary  phase  points 
associated  with  these  ray  trajectories  (i.e.,  are  obtained  from  saddle  point 
conditions  which  specify  the  ray  structure).  Very  little  numerical  d3ta  was 
obtained  due  to  limited  comput3t ional  capability  and  a  non-optimal  choice  of 


complex  Integration  paths.  Furthermore,  the  tangential  magnetic  field  on  the 
wedge  surface  and  the  far  fields  were  not  formulated.  After  elucidating  and 
extending  the  theorv  to  include  all  tangential  field  components  on  the  wedge 
surface  and  formulating  a  more  efficient  numerical  scheme  which  reduced  com¬ 
putational  time  by  one-to-two  orders  of  magnitude,  the  wedge's  surface  field 
and  directive  gain  were  show  to  possess  appropriate  physical  characteristics. 
These  include  wedge  surface  fields  whose  intensities  peak  while  phase  varia¬ 
tions  indicate  the  launching  of  radiation  and  end-fire  radiation  patterns  which 
become  more  directive  as  the  wedge  length,  relative  to  wavelength,  is  increased. 

Since  we  are  interested  in  a  dielectric  slab  waveguide  feeding  a  di¬ 
electric  wedge  antenna,  the  modal  solution  of  the  infinite  dielectric  slab  is 
relevant  and  is  reviewed  in  Section  2.  In  Section  3,  a  ray-optical  solution 
in  a  wedge  geometry  is  presented.  It  is  necessary  to  be  familiar  with  this 
solution  in  order  to  understand  the  construction,  in  Section  4,  of  the  plane 
wave  integral  representation  for  the  field  in  a  dielectric  wedge.  Numerical 
results  are  discussed  in  Section  5,  wherein  comparisons  are  made  between  the 
above  plane  wave  integral  approach,  the  local  mode  approach  of  Schering  [11], 
and  experimental  near  field  measurements  made  by  Maurer  and  Gopen  [21]. 


2.  THE  DIELECTRIC  SLAB 


Since  the  intended  purpose  of  this  study  is  to  examine  the  guidance 
and  radiation  properties  of  a  dielectric  wedge  fed  by  a  dielectric  slab 
waveguide,  it  is  convenient,  first,  to  review  the  field  structure  of  the 
infinite  dielectric  slab  in  detail. 

Exact  Modal  Solution 

The  guided  mode  structure  of  the  dielectric  slab  waveguide  is  well 
known  [22,23].  Consequently,  pertinent  results  will  be  presented  without 
detailed  derivations.  The  geometry  of  an  infinite  slab  waveguide  is  de¬ 
picted  in  Fig.  1.  Region  A  (  |x|  <  d,  jy|  <  «  )  is  occupied  by  a  non¬ 
magnetic  dielectric  of  relative  permittivity  while  regions  B~  (  |x|  >  d 
| y|  <  00  )  are  free  space.  Since  the  geometry  is  invariant  with  the  z-coor 
dinate,  the  spatial  dependency  of  the  source  free  field  structure  it 
supports  is  taken  to  be  independent  of  the  z-direction.  A  time  dependence 
of  expf-iait]  is  assumed  and  suppressed.  Of  interest  is  the  lowest  order 
even  TE  -  mode. 


For  modes  TE  to  y.  Maxwell's  field  equations  yield  the  component  set 
(Ez>  H^,  Hy) .  The  electric  field  component  satisfies  the  two-dimensional 
reduced  wave  equation 


Hu*  wavenumber  k 


-  ,e  i  (j  i  reduces  to  tlie  free  space  wavenumber  when 
o  o  r 

,  the  relative  dielectric  constant,  is  unltv;  11  ami  are  the  perme- 
r  o  o 

ability  and  permittivity  of  free  space.  Using  the  method  of  separation 
of  variables  and  considering  waves  travelling  in  the  ♦v-direet Ion  onlv, 
we  obtain  as  a  solution  to  (2.1)  the  electric  field  In  the  dielectric 
( j  x|  <  d)  to  be: 


E  "  A  cos  (k  x)  e 
?.  x 


Ik  v 
v 


k  ik  v 

V  V 

H  "  — A  cos  Ik  x)  e 
x  uin  x 

o 


12.  . 


where 


k  ik  y 

ll  ■  - — —  A  sin(k  x)  e  ' 
v  imp  x 

o 


k’  +  k'  -  k‘  “  to1:]  t  ) 
x  v  o  o  r 


l 2 . 2a' 


while  the  field  structure  in  free  space  (|x|  s  d)  can  be  shown  to  take  the 
form 


F.  •  A  cos(k  d)  e 
>•  x 


— Otx(|x|-d)  tk^v 


k  — !  x  -d)  ik  v 

v  x  1  v 

H  *  — •—  A  co8(k  d)  o  o 

X  (i)il  X 

o 

.  “01  ( |  x !  -d)  ik  V 

X  V^X  o  v 

H  "  t — r  •; -  A  cos  k  d  e  e 

v  I X I  imp ,  X 


(2.  " 


wi  th 


->'■  +  k'  -  k’  «  m'p  r 

x  V  O  O  O 


(2.  la) 


The  amplitude  constants  of  E^  In  (2.2)  and  (2.3)  have  been  chosen  to 
satisfy  continuity  of  the  E^  component  across  the  air-dielectric  inter¬ 
faces  at  x  •  ±  d.  Continuity  of  the  tangent  magnetic  field  components 
H  at  x  “  t  d  yields  the  eigenvalue  equation 


tan  k  d 
x 


(2.4) 


In  addition,  the  above  continuity  conditions  require  the  y-direeted  pro¬ 
pagation  constant  ky  in  the  dielectric  and  in  the  free  space  to  be  iden¬ 
tical. 

Bv  replacing  the  cosine  function  in  (2.2)  by  exponential  functions, 
the  field  solution  can  then  be  interpreted  as  a  superposition  of  the  two 
plane  waves,  one  progressing  upward  toward  the  interface  at  x  *■  +  d  and 
the  other  progressing  downward  toward  the  interface  at  x  *  -  d.  Hence, 
we  may  write 

E  *  E+  +  K~  . 

z  z  z 

with  (2.5) 

itk  x  +  ik  v 

E"  »  (A/2)  e  *  V 

z 

It  is  apparent  from  (2.2a)  that  these  plane  waves  impinge  on  the  inter¬ 
faces  at  angles  w^,  measured  relative  to  a  surface  normal,  given  bv 

tan  w*  *  k  /k  (2.b) 

y  x 


Since  the  modal  field  takes  the  form  of  a  superposition  of  two  plane 
waves,  one  incident  and  the  other  reflected  frimi  the  air-dielectric 


-<■>- 


One  can  show  that  the  exact  modal  solution  given  by  (2.2)  can  be  de¬ 
rived  using  the  concepts  of  ray  optics  [lb, 13, 24],  which  are  formally  pre¬ 
sented  in  Appendix  A.  For  propagating  in  the  +v  direction,  the  ray-optical 
solution  takes  the  form  of  the  two  plane  waves  of  equal  amplitude  given  bv 


eq.  (2.5).  The  ray  structure  in  the  slab  is  depicted  in  Fig.  2,  wherein 
rays  progressing  toward  the  upper  wedge  wall  (x  “  d)  are  identified  by  the 


normalized  phase  function  and  the  rays  progressing  toward  the  bottom 
wall  are  identified  by  the  phase  S0  (see  Appendix  A). 


3.  THE  DIELECTRIC  WEDGE 


While  ray  optical  and  full  wave  solutions  are  identical  for  the  dielec¬ 
tric  slab  problem,  such  agreement  is  not  attained  for  the  dielectric  wedge. 

Maurer  and  Felsen  [16]  have  effected  ray  solutions  for  the  wedge  geometries 
which  possess  constant  surface  reflection  coef f icients.  This  enabled  the  ei- 
conal  equation  of  ray  optics  (see  eq.(A4))to  be  solved  by  the  separation  of 
variable  technique,  as  was  done  for  the  dielectric  slab  [16,24].  Unfortunately, 
this  constraint  on  the  reflection  coefficient  is  not  physically  realizable  for 
the  dielectric  wedge.  In  the  wedge,  complex  diffraction  fields  of  the  same  or¬ 
der  of  magnitude  as  the  ray  optical  fields  occur.  These  diffraction  effects  are 
due  to  the  occurrence  in  the  integral  representation  of  double  saddle  points 
and  multiple  branch  points,  which  are  closely  related  to  caustics  and  multiply- 
reflected  lateral  waves.  A  caustic  is  an  envelope  of  a  system  of  real  rays. 

Lateral  waves  are  waves  that  transmit  energy  along  a  boundary  between  two  me¬ 
dia  when  a  wave  in  the  optically-denser  side  hits  the  interface  at  the  critical 
angle;  as  the  wave  travels  parallel  to  the  boundary  in  the  optically-thinner 
medium,  it  sheds  energy  back  into  the  optically-denser  material  by  refraction  [25]. 

It  has  been  recognized  that  the  simpler  approach  of  ray  optics  can  pro¬ 
vide  important,  though  incomplete,  information  about  asymptotic  fields.  Con¬ 
structs  of  ray  optics  and  results  of  a  full-wave  analysis  are  rigorously  re¬ 
lated  [27].  In  particular,  the  ray  fields  provide  arguments  of  the  functions 
that  enter  into  the  uniformly  asymptotic  descriptions  of  exact  field  solutions. 
Thus,  the  method  of  ray  optics  provides  a  convenient  starting  point  for  study¬ 
ing  the  wave  guidance  and  radiation  properties  of  the  dielectric  wedge.  Con¬ 
sequently,  a  review  of  the  ray  optical  solution  is  warranted. 
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w  with  these  ray  trajectories  such  that  cosw  =  a/r.  For  observation  points 
r  <  a,  no  real  rays  exist.  In  this  region  we  have  complex  evanescent  ray 
solutions  [lb].  By  requiring  wedge  rays  at  r  =  r^,  ,  9  =  0,  i.e.,  at  the 
interface  between  the  slab  and  the  wedge  to  evolve  continuously  from  cor¬ 
responding  slab  rays  of  incidence  angle  w  =  w^  =  w^,  (see  Fig.  3  and  eq.(2.6)), 
geometric  construction  stipulates  that 

WT  =  cos  l(a/rT>  .  (3.8) 


It  is  apparent  that  each  modal  solution  in  the  slab  corresponds  to  a 

set  of  rays  in  the  wedge  where  each  set  is  associated  with  an  eigenvalue 

a  =  a  which  defines  a  cylindrical  caustic  r  =  a  .  Since  we  are  considering 
m  m 

only  the  lowest  order  mode  in  the  slab,  we  need  only  consider  one  modal  set 
in  the  wedge  or  a  =  r^cosw^. 


To  complete  the  ray  optical  description  of  the  field  structure  in  a 
dielectric  wedge  with  surface  impedances  which  insure  constant  reflection 
coefficients,  the  amplitudes  of  the  wedge  rays  must  be  ascertained.  These 
amplitudes  were  found  by  requiring  either  conservation  of  energy  in  a  ray 
tube  [20,24]  or  by  solving  the  lowest  order  transport  equation (A5) [ 16 ] .  Con¬ 
sidering  only  rays  propagating  in  the  (-r)  direction,  i.e.,  toward  the  tip 
and  assuming  =  -C,,  =  C  in  (3.5),  the  electric  field  can  then  be  shown 
to  be  expressed  by 
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The  eigenvalues  a  -  a  are  specified  bv  (3.6)  and  are  related  to  slab  ray 
directions  by  (3.8).  The  arbitrary  constant  c  in  (3.9)  follows  directly 
by  requiring  continuity  of  phase  across  the  transition  region  at  r  -  r  , 

3  ■  0  in  the  limit  of  small  8  (see  Fig.  3). 

B .  Plane  Wave  Integral  Solutions  Whan  Ref leet ion  Coefficients  Are 

Angular  Dependent 

A  modal  ray  description  of  the  field  in  a  dielectric  wedge  with  constant 

reflection  coefficients  at  the  air-dielectric  interfaces  and  for  propagation 

towards  the  tip  in  the  region  r,^  <  r  <  a,  -  8/2  9  <  8/2  has  been  obtained 

(see  3.9).  It  is  apparent  from  (3.9)  that  the  field  amplitude  and  hence 

energy  densitv  become  infinite  as  r  -*■  a.  This  behavior  is  also  indicated 

bv  the  convergence  of  rav  tubes  at  the  cylindrical  caustics  r  »  a  .  Thus, 

m 

simple  geometric  optic  arguments  break  down  and  indicates  that  a  more  uniformly 
asymptotic  formulation  of  the  modal  field  solution  in  the  wedge  is  required. 

In  the  integral  formulation  by  Maurer  and  Siwiak  [19,20],  multiple  la¬ 
teral  wave-type  diffraction  effects  are  indicated  by  the  occurrence  of  multiple 
branch  points  in  the  angularly  dependent  reflection  coefficients.  Attempting 
to  isolate  the  individual  lateral  waves  and  their  reflections  was  not  possible 
since  the  interaction  of  the  branch  point  effects  was  found  to  be  important , 
i.e.,  they  could  not  be  treated  as  isolated  branch  points.  The  overall  effect 
could  only  be  accounted  for  by  numerical  integration  techniques.  A  reformula¬ 
tion  of  these  plane  wave  integral  representations  is  given  below. 

To  justify  the  postulated  form  of  the  integral  representations  for 
the  modal  field  on  the  surface  of  a  dielectric  wedge  (see  (3.1**^  and  i3.35)), 
consider  two  plane  wave  solutions  to  the  two-dimensional  scalar  wave  equa¬ 
tion  which  are  assumed  to  propagate  in  the  +y  direction. 


+  -  ,  +  jkxx  -  -jkx-x  jkv-v 

E=E+E=(Ae+Ae  *  )  e  y 
z  z  z 


(3. 10) 


where  k  *  /  k‘  -  k'  >  0.  By  transforming  this  solution  into  polar  coor¬ 


dinates  via  the  relations 


x  =  r  cos  f 


y  =  r  sinC 


k  =  k  cos  a 
x 


k  =  k  sin  a 
v 


we  obtain 


E  =  E+  +  E"  =  A+  e‘jkr  C°S  (n+a_0)  +  A~  e'jkr  co8<***>  (3.11) 

z  z  z 

Eq.  (3.11)  represents  two  plane  waves  which  progress  in  the  directions 

* 

$  =  n  +  a  and  *>  «*  -a  toward  the  origin  in  xy-space,  respectively. 


A  more  general  solution  to  the  wave  equation  can  be  constructed  by 
superposition  f 26 J .  Hence,  we  may  write 


e+-  r 

z  j 

ai . 

/  s„ 


\  -jkr  cos  (H+a-^i)  , 

A  (a)  e  J  r  da 


A~  (a)  e"jkr  cos^  da 


(3.12a) 


(3.12b) 


with  A~ (a)  an  analytic  function  dependent  on  a.  E^  represents  a  bundle  of 
plane  waves  propagating  in  the  directions  a^  £  a  £  toward  the  upper 
wedge  air-dielectric  interface;  E^  identifies  a  set  of  plane  waves  progres¬ 
sing  toward  the  lower  wedge  wall  in  the  directions  a.-,  £  a  <_  So  • 


Specification  of  wave  motion  toward  the  origin  is  relevant  in  polar  coordi¬ 
nates  since  our  main  concern  is  with  the  effects  of  wave  guidance  toward  the 
tip  of  the  wedge. 
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In  an  unbounded  domain,  solutions  to  the  wave  equation  in  polar  coor¬ 


dinates  are  of  the  form  f  26  ] 

E  =  Z  (kr)  elv°  ,  (3.13) 

Z  V 

where  v  is  a  separation  constant  and  Z  (kr)  represents  a  general  cylinder 
function.  For  waves  approaching  the  origin,  Z  becomes  a  Hankel  function 
of  the  second  kind  whose  integral  representation  in  the  complex  w' -space 
is 


-5 


r 


+ikr  cos  w'  ika(w'  -  11/2)  ,  , 
i  e  dw' 


(3.14) 


where  v  has  been  set  equal  to  ka  and  the  path  P,  is  defined  in  Fig.  5. 


The  ray  optical  interpretation  of  the  asymptotic  form  of  solution 
(3.13)  with  given  by  (3.14)  and  v-positive  is  depicted  in  Fig.  6.  As 
illustrated,  rays  3re  trajectories  orthogonal  to  the  wave  fronts  and  tangent 
to  the  circular  caustic  r  =  a.  Since  our  wave  bundles  (3.12)  also  approach 
the  origin  and  we  expect  product  separability  as  in  (3.13),  let  us  assume 
that  the  angular  dependent  coefficients  in  can  be  expressed  as 


A  (oi)  *  B  ejka(a  H/-)  ^  g-  independent  of  at. 
Substitution  of  (3.15)  into  (3.12a)  and  (3.12b)  gives 
s-Bl 

_+  _  „+  -jkr  cos(H+a-iJ>)  +  jka(ct-II/2) 

t  "  w  0  U  Ot 

2  J 

“*  ?, 

E*  -  9'  j  '  e"Jkr  c0*<*rt+>  +  Jkata-a/2)  J0, 

2  °2 


(3.15) 


(3. 16a) 


(3. 16b) 
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From  Fig.  3,  it  is  clear  that  the  coordinate  ,s  is  more  convenient  than  the 
?  coordinate.  Setting  8  ■  $  -  3a/ 2  reduces  (3.16a)  and  (3. lob)  to 


+  + 

E  -  B 


^-jkr  sin(i-'^)  +  Jka(<i-1/2)  j 


(3. 17a) 


E  -  B 


-jkr  sin(a+8)  +  jka(a-H/2) 


(3.17b) 


As  in  (3.13),  the  '-dependence  can  be  isolated  by  introducing  the  transfor¬ 
mations  w  »  a  -  8  in  E+  and  w  »  a  +  3  in  E  ;  hence 

2  z 


3  —3 

.i  .rr/1  /  l  -jkr  sinvf  jkaw, 

_+  +  -jka  (11/2-8)  /  e  dw 

E  -  B  e  / 

=  > 


(3.13a) 


-jka(!I/2+8)  ^"Jkr  SlnWf 

E^  -  B  e  J 

a -,+8 


(3. ISb) 


At  the  upper  wedge  wall,  i.e.,  when  8  -  8/2,  the  ray  bundle  E^  in  (3.13) 
reflects  into  the  ray  bundle  E^  as  shown  in  Fig.  7.  From  the  geometry  we 
see  that  ji  varies  in  (3.18a)  between  *  w^,  -  8/2  and  8.,  -  w^.  +  8/2,  where 

Wj  is  equal  to  the  angle  of  incidence  of  the  plane  waves  in  the  slab,  and 

33  ,  3 

in  (3.13b)  between  -  —  and  -  —  .  Hence  over  region  1,  (see 

Fig.  7  for  specification  of  region  1  on  the  upper  wedge  wall),  the  field  is 
given  by 


Ezl  ■  Ezl  +  Ezl 


(3. 19) 


where 
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E+.  -  B+  ' 

21  vs 


ikaw  -  ikr  sinw 
e  dw 


(3. 19a) 


E-  .  B-  e-jka(n+S)/2  J 

21  e-s 

T 


ikaw  -  ikr  sinw 
e  dw 


(3. 19b) 


Furthermore,  B  =  B  .  This  follows  from  symmetry  by  noting  that,  for  ex¬ 
citation  of  the  wedge  by  a  symmetric  slab  mode,  the  field  over  the  top  and 
bottom  surfaces  of  the  wedge  (9  =  +3/2  and  -3/-,  respectively)  must  be  iden¬ 


tical.  Consequently,  the  incident  wave  bundle  E  ^  illuminating  region  1  on 
the  upper  wedge  surface  equals  the  field  which  illuminates  region  1  on 
the  lower  wedge  surface,  and  both  are  superpositions  of  plane  waves  which 
span  the  angular  range  w^  -  g/2  <  w  <  w +  g/2.  Hence,  from  eqs.  (3.18a) 
and  (3.18b),  B+  =  B  =  B^  and  we  may  write  for  the  field  over  region  1  on 
the  upper  wedge  wall  (9  =  8/2): 


E  ,  =  B  e 
zl  l 


-jka(II-S)/; 


[r1  <■  *  * 


JkW<w>  dw 


1  9_-8 
^  T 


(3.20) 


where  the  reflection  coefficient  on  the  upper  wedge  wall  is  given  by 


r+  ,  E~  /  E+  -  ~jka8  =  r 
1 1  zl  '  zl  6  “  ‘t 


(3.21a) 


W(w)  e  aw  -  r  sinw 


(3.21b) 


Eq.  (3.21a)  also  follows  from  the  modal  resonance  relation  (3.6)  for  the 
lowest  order  mode  (ra  =  0).  Furthermore,  for  small  apex  angles  6  it  follows 
from  (3.8)  and  the  relationship  sin  3/2  ~  8/2  =  d/rT  that  expression  (3.21a) 


reduces  to  eq.  (.2.7)  with  F~  given  by  the  Fresnel  polarization  dependent 
reflection  coefficient  (2.9). 


The  integral  representations  for  the  magnetic  field  components  follow 
directly  from  Maxwell's  source  free  equations 


1  He,  .  h. 


l 


-fr  E. 


r  j  r  '  9  z  b  jwii  jr  z 

•'ll  ^  n 


(3.22) 


Applying  (3.22)  to  the  integrals  (3.17a)  and  (3.17b)  and  following  the  same 
procedure  which  led  to  obtaining  (3.18a)  and  (3.18b)  results  in  giving  for 
the  magnetic  field  components 


H  -  H+  +  H"  ,  H,  =  H+  +  H” 

r  r  r  9  9  9 


(3.23) 


with 


,-6,-9 


jkW(w)  .  -j  ka  (11/ 2—0 ) 
cosw  eJ  dw  e  J 


c^-0 


(3.23a) 


H  =  i- 


and 


,  30+9 


jkW(w)  J  -jka(H/2+9) 
cosw  eJ  dw  e 


a^+6 


/  +  y6rd 

I*-  '  sinw  ejkW<w)  dw  e-jka(II/::-e) 

V  n 

■v9 


(3.23b) 


(3.23c) 


9  1  n 


-S2+0 


jkW(w)  ,  i  -jka(H/2+6) 
sinw  eJ  dw  1  e  J 


,-K? 


(3.23d) 


where  W(w)  is  given  by  (3.21b).  Hence,  for  9  =  3/2  and  illumination  over 
region  1,  (3.23)  reduces  to 
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the  wedge  Coward  the  apex.  Consider  the  situation  of  once  reflected  beams 

and  illumination  of  region  2  on  the  wedge  air-dielectric  interface  (see  Fig. 

7).  From  geometric  construction,  upgoing  plane  waves  which  illuminate 

region  2  are  confined  to  the  angular  directions  w^  -  36/2  <  a  <  -  6/2, 

whereas  reflected  waves  E  .  are  confined  to  the  directions  w_  -  56/2  <  a 

zz  L 

<  Wj,  -  36/2.  As  a  consequence,  eqs.  (3.18a)  and  (3.18b)  take  the  form 


E+,  =  fit  e-jka(II"3)/2  f  ejkW(w)  dw 
z2  2  j 

w  -26 


(3.26a) 


E",  -  Bt  fT  '  dw 

z2  2  J 

wT  -26 


(3.26b) 


where  W(w)  is  given  by  (3.21b).  Hence,  at  9  =  8/2  along  region  2  in  Fig. 
7,  the  surface  illumination  from  the  beam  progression  toward  the  tip  is 
given  by 


w  —  8 

_  _+  .  P—  n  -jka ( IT—  6 ) / 2  <  T  ,  r+.  jkW(w)  , 

E9=E-+E0=B.  eJ  I  (l+i„)eJ  aw 

zZ  zZ  zZ  Z  j  l 

wt-28 


(3.27) 


with  reflection  coefficient  T 2  given  by 


r+  =  E~  /E+  =  e~^ ka3  = 

*2  -  “z2'  z2  ® 


(3.27a) 


and  symmetry  requiring  that  =  B?  =  B?. 


In  order  to  relate  B«  to  B, ,  the  observation  is  made  that  E  ,  waves 

2  1  zl 

incident  on  region  1  of  the  lower  wedge  air-dielectric  interface  span  the 


angular  directions  w^,  -  6/2  <  a  <  w^  +  8/2  and  reflect  as  E^  in  direction 
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E*1  -  »1 


1 


| kW 


i  K/S.O 


t :sm 


Vho  reflection  coefllclont  along  region  1  ot  the  lower  .» I  r-»it  o  lee  t  r  i c  ln- 
tertace  Is  now  ileflneii  .is 


r;  ’  ",  - 


( i.  :»> 


who  r  o  I*  has  boon  Int  rodviceii  in  fl..'lal.  On  replacing  IV,  in  v'..’bhl  bv 
using  13.2*0  ami  combining  K*  ,  witti  E.  ,  ono  obtains  tor  t  ho  surface  field 

2  •-  2  i 

along  rogiv'n  on  the  upper  wedge  wait  t  ho  expression 


K,>  «  »,  * 


-)ka(ll-t')/:  r  ’  l"  f!  1 '  JkWiwI 


t  l.  101 


w  -  ’  g 

r 


Again  allowing  the  reflection  coefficients  to  be  angular  itoponilont  ami  ob¬ 
serving  that  it  rj  *  !'iw^  then  !' ^  must  equal  1 ' t w  V  tO ,  eq.  i  h"  becomes 

E  ,  -  8.  f  1'A"'„+  ii)  ,,  +  r (w )  ,  o,kVKw'  ,lw  n.m 

**■  y  r 

wr-.b 

where  Piw^  Is  given  by  y  J . ,J  A  with  v>'^  replaced  bv  w. 


8v  continuing  the  above  process  of  tracking  beam  reflections  from 


the  wedge  walls,  it  evolves  that  the  surface  field  along  the  n 


region 


lit  the  upper  wedge  w.i  1  1  t'.'in  be  approx  l  mat  ed  I'v  t  ho  iiitonf.il  express  ion 


l )  6 

„  -lkaUl-bl/2  /  n-1  !' (w  (•  mo')  IklJtwl 

K,n  “  "l  °  J  [  1+1’ (w)  ]  11  ,.n- 1  Jw 

wl-nb  »<“  1  '  r 


i  i.  121 


Thus,  along  the  plane  0  -  d/2  between  the  apex  r  -  0  anil  the  transition  ft 


Hli’ii  r  »  r.j,,  the  surface  electric  field  Is  given  by 

N 

k  -  e  k  .  ii.m 

z  zn 

n-  1 


where  K  Is  spool  fled  bv  (\.  121  and  w_  -  Nd.  The  Integer  N  Identities 
zn  1 

the  number  ol  regions  along  the  upper  wedge  wall  that  are  dlrootlv  illuminated 

by  the  plane  waves  as  they  progress  toward  the  tip  and  are  multinlv  refleoted 

from  the  wedge  alr-dle  lectrlo  Interface.  Note  that  when  n  -  N  w  -  w.^.  - 

nd  •  0;  hence,  the  final  wave  constituent  In  F  is  normally  incident  on  t ho 

zn 

upper  wedge  wall.  As  was  pointed  out  In  l'a|,  this  plane  wave  is  tangent 
to  the  caustic  at  r  ■  a,  0  -  d/2. 


For  n  N  N,  the  Integral  (.1.  121  describes  surface  illumination  for  waves 
progressing  back  toward  the  slab  region.  When  n  m  2N ,  w  -  -w  ;  the  asso¬ 
ciated  plane  wave  constituent  Illuminates  region  l  on  the  upper  wedge  wall 
at  observation  point  ir^.,  -11/21  and  Is  incident  at  an  angle  of  -w^,  measured 
to  the  right  ot  a  unit  normal  to  the  wedge  wall  |  2-t  ] .  Hence,  the  surtace 
Illumination  considering  both  right-going  and  left-going  waves  is  given  bv 
the  expression 

2N 

F.  -  E  F  .  0„  -  Nfl  ,  l  i.  Is l 

z  zn  1 

n»  l 

which  results  in  Integers  being  evaluated  numerical lv  between  the  limits 
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-  w  <  w  <  w  _ 
T  T 


The  remaining  field  components  can  be  inferred  from  (3.25)  and  (3.32). 
Hence, 


2N 

1  =  Z  H 

r  rn 

n=  1 


2N 
Z  H 
n=l 


9n 


(3.35) 


with 


n-1 


,  _l  e-jka<n-6)/2  f  £os  (1  .  r(v))  Y  d„ 

rn  n  ;  o-i  r”'1 

1  (3.35a) 


-ng 


and 


wT-(n-l)B 

n  -  _i  -jka(II-B)/2  f 

H0n  n  e  J 

w„-nB 


.  r,  „ ,  x  ,  r(v+mB)  ikW(w) 

sin  w  [1  =  i  (w)  ]  II - ; -  eJ  dw 

m-L  r"-1 

T  (3.35b) 


The  above  expressions,  eqs.  (3.34)  and  (3.35),  are  assumed  to  give  the 
surface  field  at  observation  points  r<a,  i.e.,  in  the  region  shielded  by 
the  caustic.  The  justification  for  this  assumption  is  that  complex  dif¬ 
fraction  effects  take  place  when  the  incidence  angle  w  passes  the  critical 
angle  9^,  which  is  given  by  0c  =  sin  ^(1/  i/£TJ)  and  corresponds  to  a  branch 
point.  The  integration  paths  of  (3.34)  and  (3.35)  encounter  multiple  branch 
points  that  are  introduced  by  the  multiple  reflection  coefficients  in  the 
integrand.  These  branch  points  establish  lateral  wave-types  that  constitute 
the  dominant  field  within  the  caustic  [24]. 


Although  the  integration  paths  lie  in  the  range  -w  <  w  <  w^,  numerical 
evaluation  of  (3.34)  and  (3.35)  as  well  as  experimental  evidence  have  in¬ 
dicated  that  no  measurable  energy  is  reflected  back  into  the  slab.  This 
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THE  FAR  FIELD 


Using  the  theory  presented  in  Section  3,  the  tangential  electromag¬ 
netic  field  components  can  be  calculated  at  a  suitable  number  of  points 
along  the  surface  of  a  given  wedge.  From  the  results  of  Section  2,  the 
evanescent  field  external  to  the  slab  can  be  determined  along  a  hypothetical 
wedge  boundary  extending  out  to  infinity  radially,  as  diagramatically  il¬ 
lustrated  in  Fig.  8.  These  radial  lines,  in  conjunction  with  a  circular 
path  at  infinity,  enclose  a  domain  D,  occupied  by  free  space  which  ic 
external  to  the  dielectric  slab-wedge  radiator.  Since  domain  D  contains 
no  sources,  Ki  rchhof f-Huygen  diffraction  integral  can  be  applied  to  find 
the  fields  radiated  into  the  region  D  in  terms  of  the  tangential  fields 
on  tlie  boundary  of  D,  designated  3D^,  i  =  1,2  and  in  Fig.  8. 

In  free  space,  the  radiated  electromagnetic  field,  E^(r)  satisfies 
the  two-dimensional  wave  equation: 

(V2  +  k2)  E  (r)  -  0  (4.1) 

o  z  — 

where  Ez(r.)  is  the  radiated  field,  r  is  the  position  vector  to  some  point 
P  in  D  where  the  far  field  is  to  be  determined,  and  k^  is  the  free  space 
wavenumber.  Let  G^  satisfy  the  equation 

(V2  +  k2)  G  (r,  r')  -  -  5(r  -  r')  (4.2) 

o  o  —  —  —  — 

where  &(r_  -  _r'),  the  Dirac  Delta  function,  represents  an  inf initessimal 
point  source  of  unity  magnitude  at  r  ■  r ' .  By  first  multiplying  (4.1) 
by  and  (4.2)  by  -E^(r'),  subtracting  the  two,  then  integrating  the  re¬ 
sultant  expression  over  the  surface  enclosing  D  and  using  Green's  second 


Identity  in  two  dimensions,  it  evolves  that  the  far  field  can  be  obtained 


between  the  normal  n  to  the  wedge  surface  and  position  vector  R.  The 
normal  derivative  of  the  Green's  function  can  then  be  shown  to  take  the 
form 


3G  ,  3G 

_ o  _ o  _ 

3n  r  3<t 


k  cos(u)G 
o  o 


(4.7) 


with  Go  given  by  the  asymptotic  form  (4.5). 

From  Maxwell's  equation  in  cylindrical  coordinates,  the  normal  deriva¬ 
tive  of  the  tangential  electric  field  is  related  to  the  tangential  magnetic 
field  by  the  expression 


3E  ^  .  3E 

2  1  _ z 

-r —  -  n  .  VE  =  —  -r— 

3n  —  z  r  30 


-July  H 
or 


(4.8) 


Combining  (4.7)  and  (4.8)  with  (4.3)  gives  the  radiated  far  electric  field 


E 

z 


(r) 


[cos(u)E2  -  no  H^] 


dr' 


(4.9) 


In  (4.9),  the  integration  ranges  over  the  two  segments  3D^  and  3D,  defined 
in  Fig.  8,  Gq  is  given  by  (4.5)  and  (kQ,  r)Q)  are  the  wavenumber  and  the  in¬ 
trinsic  wave  impedance  of  free  space.  The  integration  of  (4.9)  is  carried 
out  along  both  wedge  surfaces.  As  will  later  be  discussed  in  greater  de¬ 
tail,  Romberg’s  method  was  employed  in  the  analysis  of  actual  wedge  prob¬ 
lems  . 
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s.  NUMERICAL  RESULTS 


The  wedge  antenna  is  completely  characterised  by  three  parameters: 
kod.  the  relative  half-width  of  the  slab  with  respect  to  the  free  space 

wavenumber  k  ,  k  ?.  .  the  relative  length  of  the  wedge  and  »'  .  the  relative 

o  o  v  r 

permittivity  of  the  antenna  material.  These  values  are  varied  so  as  to 

determine  the  effect  each  has  on  both  the  surface  field  and  the  far  field. 

In  addition,  surface  and  far  fields  are  calculated  for  two  antennas  with  a 

relative  permittivity  of  silicon,  »  12.  For  these  two  antennas,  our 

theoretical  results  are  compared  with  an  alternative  approach  based  on  the 

local  mode  theory  [11].  The  silicon  material  with  t'  »  12  is  currently  un- 

r 

der  study  for  use  in  millimeter-wave  transmitters  and  receivers  being  de¬ 
veloped  by  the  U.S.  Army.  The  antenna  parameters  examined  are  listed  in 
Table  1,  with  identifying  letters  A  through  I,.  Also  included  in  this  table 

are  the  values  of  w_,  3  and  k  a.  The  angle  w_  is  the  incident  angle  for 
i  o  “■  i 

the  plane  waves  in  the  slab  which  excite  the  wedge.  Recall  that  the  range 
-wt  £  w  <  w^  specifies  the  integration  limits  for  the  wedge  surface  fields. 
The  parameter  b  is  the  apex  angle  of  the  wedge  and  k^a  is  the  location  of  the 
caustic  relative  to  the  free  space  wavenumber.  Antennas  J,  K,  and  L  are  of 
particular  interest  because  of  the  availability  of  empirical  data  collected 
by  Maurer  and  Gopen  [21].  The  material  used  in  the  experiment  was  rexolite 
with  a  relative  permittivity  f  -  2.5h. 

In  Section  4  and  in  eq.  it  was  shown  that  a  far  field  evaluation 

requires  first  calculating  the  field  components  E,  and  H  along  the  planar 
surfaces  30^  and  ?D,,  identified  in  Fig.  3.  Physically,  these  surface  field 
components  must  be  continuous  along  and  IP,.  Hence,  we  imposed  the  con- 


dition  chat  the  electric  field  on  the  surface  of  the  dielectric  wedge  be  equal 
to  the  electric  field  in  the  evanescent  region  of  the  slab  waveguide  at  the 
trnasition  point  r  »  r^,  0  =  8/2,  i.e.,  we  equated  (3.33)  to  (2.3).  The 
small  discontinuity  in  H  across  the  transition  point  was  found  to  be  accep¬ 
table  (see,  for  example.  Fig.  13  through  Fig.  16).  In  addition,  it  was  con¬ 
venient  to  normalize  the  surface  wedge  field  by  setting  the  multiplicative 
factor  in  (3.32)  Bj  exp  [  —  j ka (IT  -  3/2)]^  to  unity. 

A.  The  Surface  Field 

The  relative  magnitude  (in  dB)  of  the  surface  electric  and  magnetic 
tangent  fields  as  well  as  their  respective  phases  are  plotted  versus  the  nor¬ 
malized  coordinate  k^y'  in  Figs.  9-16.  Values  are  plotted  from  the  tip 

(k  y'  =  k  r  ),  past  the  transition  region  (k  y'  =  0)  and  extend  into  the 
o  o  r  o 

evanescent  region  of  the  slab  waveguide.  The  phase  difference  at  the  apex 

of  each  antenna  is  normalized  to  read  411  radians.  In  Figs.  9a  and  9b  and 

Figs.  13a  and  13b,  the  wedge  length  is  increased  while  k^d  and  Gr  are  held 

fixed.  Figs.  10a  and  10b,  together  with  Figs.  14a  and  14b  show  the  effects 

of  changing  the  relative  slab  width  k  d  while  Figs.  11a  and  lib,  with  Figs. 

o 

15a  and  15b  depict  the  surface  field  as  G^  is  changed.  In  Figs.  16a  and 
16b,  the  surface  field  for  two  antennas  with  Gf  =  12  is  presented.  How 
these  changes  affect  the  far  field  will  be  discussed  shortly. 


Of  particular  interest  is  the  phase  velocity  v  ^  ,  along  the  wedge 

surface  in  the  y’-direction  (forward  direction).  In  order  for  a  wave  to 
be  launched  from  the  wedge  in  this  direction,  the  corresponding  component 
of  the  wavevector  k^,  must  equal  the  free  space  wavenumber  k^  (this  follows 
from  Snell's  law).  The  wavevector  k  is  defined  as  the  spatial  gradient  of 
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the  phase  of  the  wave  function  (in  our  case,  E^).  Since  k^,  »  v'  .  k, 
where  v'  is  a  unit  vector  in  the  positive  v' -direction,  we  may  formulate 
the  relation 

k  ,  -  k  -  -r—  (ARG(E  )  -ky'l  (5.4) 

y  o  3y  z  o' 

where  y*  *  r^  -  r'.  When  this  quantity  is  zero,  the  surface  wave  is  launched. 
Graphically,  this  is  determined  by  plotting  the  function  ArglE^)  -  k^y' 
and  noting  that  when  the  slope  is  zero  the  surface  wave  enters  the  'fast 
wave'  region  and  radiation  becomes  possible.  This  condition  is  illustrated 
in  Figs.  9b,  10b,  lib,  and  12b. 

B.  The  Far  Field 

As  the  length  k  1  of  the  antenn3  is  increased  (see  Figs.  9  and  12) 
o 

there  results  an  expansion  of  the  spatial  interval  on  the  surface  of  the 
wedge  over  which  the  surface  field  has  a  phase  velocity  in  the  y'-direction 
equal  to  the  speed  of  light.  Consequently ,  one  expects  enhanced  radiation 
capability.  In  Figs.  17  through  20,  the  directive  gain  g^  is  plotted  versus 
the  angle  which  Is  defined  in  Fig.  8.  For  two-dimensional  geometries, 

the  directive  gain  is  given  by 

gd  -  2trrSr/Prad  .  (5.5a) 

where  is  the  real  Povnting  vector  in  the  radial  direction,  expressed  in 

terms  of  the  far  electric  field  bv  S  «■  !  E  ( r )  ! 2  / n  ,  r)  »  (u  It  )  *  the 

"  r  z  —  o  o  o  o 

intrinsic  Impedance  of  free  space,  and  P^d  is  the  power  radiated  by  the 
antenna.  From  Fig.  17,  it  appears  that  a  larger  value  of  leads  to  more 


power  in  the  forward  direction. 


tics 


Bv  decreasing  the  width  k  d  of  the  antenna,  the  endfire  characteris- 

o 

are  also  improved.  Thus,  in  Fig.  13,  antenna  E  has  more  directivity 
than  either  antenna  C  or  D.  This  is  most  likely  due  to  the  fact  that  an¬ 
tenna  E  possesses  a  larger  'fast  wave'  region  along  the  wedge  surface  than 
either  antenna  C  or  D,  as  is  evident  from  Fig.  10b. 

In  Fig.  19,  the  relative  permittivity  c  is  varied.  Decreasing  this 
parameter  also  increases  the  sice  of  the  fast  wave  region  (see  Fig.  lib) 
and  thus  antenna  G  has  better  propagating  characterist ics  than  either  E  or 
F . 

Fig.  20  demonstrates  that  directive  antennas  could  be  designed  with 
silicon  (e  =  12):  however,  they  would  have  to  be  very  long  or  very  narrow. 

It  was  also  observed  that  slightly  changing  the  frequency,  while  not 

affecting  the  main  beam,  can  increase  or  decrease  the  side  and  back  lobes 

significantly,  indicating  that  a  more  desirable  combination  of  k  d,  k  1 , 

o  o 

and  can  be  found  for  a  given  application  at  a  given  frequency. 

The  directive  gain  for  antennas  H  and  I,  with  f  =  12^ are  compared 
to  patterns  generated  by  a  local  mode  approach  [11  ]  in  Fig.  20.  In  the 
end-fire  direction,  the  two  methods  give  comparable  values  for  the  gain 
function,  \Tien  q  N  30°,  the  local  mode  theory  yields  results  which  de¬ 
crease  more  rapidly  than  the  plane  waves  approach:  however,  for  these  ob¬ 
servation  angles,  the  gain  is  already  20  dB  below  its  maximum  value. 

In  Figs.  21,  22,  and  23,  experimental  results  [21]  are  compared 
with  theoretical  ones.  The  normalized  coordinate  k  v'  has  its  origin  in 

O' 

the  wedge  surface  at  the  transition  region  r'  •  rT  (see  Fig.  3)  and  is 
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directed  toward  the  wedge  tip.  The  parameter  k  is  the  free  space  wavenumber. 

o 

The  curves  are  normalized  so  that  the  experimental  and  theoretical  values  for 
the  peak  fields  coincide.  Near  field  experimental  values  were  measured  ap¬ 
proximately  one  millimeter  away  from  the  wedge  surface.  The  small  fluctua¬ 
tions  or  ripples  on  the  experimentally  determined  curves  in  Figs.  21-23  are 
most  likely  due  to  the  presence  of  standing  waves  between  the  probe  and  the 
absorbing  material  surrounding  the  wedge  and  probe  in  the  parallel-plate  ex¬ 
perimental  setup  [6].  As  is  evident  from  each  of  the  figures,  the  experi¬ 
mental  values  dip  downward  across  the  transition  region,  while  the  theoretical 
values  tend  to  tilt  upward.  There  appears  to  be  a  discrepancy  of  about  four 
to  five  dB  between  these  two  results.  It  ought  to  be  pointed  out  that  the 
location  of  the  wedge,  relative  to  the  experimentally-determined  surface 
field  values,  was  difficult  to  ascertain  from  the  experimental  setup*  and 
that  an  average  value  was  used  for  the  relative  permittivity  of  rexolite. 

The  theoretical  curves,  nonetheless,  do  exhibit  the  correct  shape  for  the 
surface  electric  field,  i.e.,  a  monotonicallv-increasing  amplitude  which 
peaks  and  then  falls  as  the  tip  is  approached. 


♦Unfortunately ,  the  experimental  setup  has  been  dismantled,  which  prevented 
duplicating  the  measurements  to  further  substantiate  their  validity. 
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CONCLUSION 


b. 

A  theory  has  evolved  which  appears  to  yield  an  approximate  solution 
for  the  surface  field  on  a  dielectric  wedge.  The  surface  field  exhibits 
meaningful  physical  characterist ics  such  as  magnitude  which  peaks  and  a 
phase  variation  which  predicts  launching.  The  analysis  is  very  general 
and  can  be  adapted  with  relative  ease  to  a  variety  of  two-dimensional 
tapers  such  as  wedges  with  curved  boundaries.  In  addition,  certain  three- 
dimensional  geometries  such  as  cones  or  pvramidal-type  structures  might 
be  amenable  to  such  analysis. 
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APPENDIX  A.  RAY  OPTICS 


The  ray  optical  solution  to  the  reduced  wave  equation 

(V:  +  k2)  u(r)  =  0  (Al) 

is  sought.  The  quantity  k  =  u)  »  u  e  e  denotes  the  wavenumber  of  a 

o  o  r 

homogeneous  lossless  dielectric  of  relative  dielectric  constant  er-  A 
time  dependence  of  exp  (-iiot)  is  assumed  and  suppressed.  To  obtain  a 
unique  solution  of  (Al).  either  the  tangential  electric  and  magnetic 
fields  at  the  air-dielectric  surfaces  must  be  continuous  or  the  tangent 
electric  field  on  the  boundary  must  be  linearly  related  to  the  tangent 
magnetic  field  (the  so-called  "impedance  boundary  condition")  which  can 
be  expressed  by  the  relation 

Z  (r)  3u(r) 

-tk  u(r)  4  _ =0  on  B  (A2) 

n  3v 


where  v  is  the  outward  normal  to  the  boundary  surface  B,  n  is  the  intrinsic 
Impedance  of  the  dielectric,  Zg(r)  is  a  surface  impedance,  and  u(r)  repre¬ 
sents  the  transverse  electric  field  component  E^. 

Following  the  development  by  Maurer  and  Felsen[lb],  a  geometric  optics 
solution  for  large  k  of  the  form 


N 

u(r)~  Z 

p-1 


ikS  (r) 

A  (r)  e  P 
P  ~ 


(All 


is  assumed.  Each  species  (denoted 
normalized  phase  S^(r).  Inserting 


by  p)  has 
(A3)  into 


an  amplitude  A^ (r)  and  a 
(Al)  and  equating  to  zero. 


the  coefficients  of  the  k2  and  of  the  k  terms  give,  respectively. 


(  VS  (r)  )2  =  1  (A4) 

P  ~ 

and 

2  V  S  (r)  .  V  A  (r)  +  A  (r)  V2  S  (r)  =  0  (A5) 

p-  p-  p-  p  — 

for  any  species  p.  Eqs.  (A4)  and  (A5)  are,  respectively,  the  eiconal  and 
transport  equations  of  geometrical  optics.  Their  solutions  provide  ampli¬ 
tude  and  phase  variations  of  the  rays.  It  should  be  emphasized  that  ray 
solutions  are  only  asymptotic  expressions  and  are  not  full  wave  solutions. 
One  way  of  explaining  the  accuracy  of  a  ray  solution  is  to  compare  it  with 
the  asymptotic  form  of  a  full  wave  solution. 


Substituting  the  ray  solution  (A5)  into  the  boundary  condition  (A2) 
yields 


Eq.  (A6) 
so  that 


ik  S  (r)  9S  (ry 

l  e  p  A  (r)  [-  1  +  Zg(r)  — 2 - 1  »  0  (A6) 

P*1 

can  be  satisfied  by  postulating  the  pairwise  vanishing  of  terms 


ik  S 

3S 

ik  S 

3S 

e  P  [Vi 

Vr>  if  -  '] 

+  e  ’  fVI 

Z  (r)  -r— ^  -  1 
s  —  3v 

^  > 

on  B 


Equating  phases  of  the  exponential  terms  in  (A7)  yields 


k  S  =  k  S  +  2Hin  on  B, 

P  q 


0 


(A7) 


(A8) 
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where  m  id  an  integer  such  that  initially  is  of  order  1.  From  the 
eiconal  equation  (A4)  and  boundary  condition  (A8) ,  it  can  be  shown  that 


3S  3S 

_ _ H 

3v  3v 


on  B 


(A9) 


The  minus  sign  enters  into  the  above  equation  because  reflection  takes 
place  at  the  boundary.  Using  (A8)  and  (A9)  in  (A7)  gives 


A  -  r  A 
9  P 


on  B 


(AlOa) 


where 


3S 


Z.<*>  3^  -  1 

3S~ 


(AlOb) 


sT  +  1 


r  is  recognized  as  a  reflection  coefficient  and  Z  *  Z  /r). 

s  s 


If  the  reflection  coefficient  is  constant  on  the  boundary  B,  it  is 
possible  to  construct  an  alternative  formation  without  destroying  the 
specular  reflection  condition  (A9).  In  this  formulation,  the  phase  of  F 
is  incorporated  into  the  phase  function  S  rather  than  being  totally  as¬ 
sociated  with  amplitude  terms  as  was  done  in  (A10).  Rewriting  (A7)  in 
the  form 

ik  S  ik  S  iargr  +  ik  S 

Aq  e  q  -  Ap  F  e  P  -  Ap|r|  e  P  ,  (All) 


with  f  defined  by  (AlOb),  leads  naturally  to  the  relations 
Aq  *  Ap|r|,  k  Sq  •  k  Sp  +  arg  F  -  2mn. 


(A12) 
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appendix  b.  computer  program  documentation 


A.  Surface  Field 

Two  computer  programs  are  used  to  compute  the  surface  field  for  the 
wedge  antenna.  Since  only  the  tangential  components  of  the  electromagnetic 
field  are  required  to  calculate  the  far  field,  only  (3.341  and  (3.35)  are 
used.  They  are  repeated  here  for  convenience. 


N 

E  =  -  E  {B  e 

z  1 

n=  1 
n-1 


.w^-ng 


■jka(n 


—  6) / 2  j  e"jkr  Sln  W  +  J  k3W  (1  +  T(w)  ), 


wT  -(n-l)S 


„  ,T(w  +  mB)  N  ,  i 
n  ( - — )  dw  } 


„n-i 


(Bl) 


m»  1 


N 


B, 


.  -  _L  -jka(IT-S)/ 

n 

n=  1 


n- 1 


■/ 


w  -n6 

e-1kr  sin  w  +  j  kaw  (1  _  f(w) 


wT  -(n-l)S 


_  ,r(w  +  mB) ,  ,  i 

II  ( - : — )  cos  w  dw  } 

.n- 1 


(B2) 


Program  'One'  evaluates  the  r-independent  part  of  these  two  integrands, 
namely: 


(l  +  T(w)  ) 


n-1 

r(w  +  mB)  -}kr  sin  w 
R  (  .n-1  e 

m»l  1 T 


1B3) 


and 


n- 1 

(i  -  r(w)  )  n 

m*>  1 


F(w  +  rat?) 
„n- 1 
1  T 


-jkr  sin  w 

£  COS  w 


( B4 ) 


as  vector  quantities  and  stores  these  vectors  in  a  disc  file  as  input 
data  to  program  'Two'  to  be  integrated  with  the  r-dependent  parts  of  the 
program.  This  method  is  employed  since  calculating  the  product  of  reflec¬ 
tion  coefficients  is  very  time-consuming  and  can  be  done  once  for  each  set 
of  antenna  parameters. 

'One'  requires  as  input:  k  d  the  relative  half  width  of  the  slab, 

o 

k^i  the  relative  length  of  the  wedge,  and  e  the  relative  dielectric  con¬ 
stant.  Five  values  are  calculated  per  6  radians  on  the  real  w-axis,  start¬ 
ing  with  and  ending  when  the  magnitude  of  the  product  term  is  less  than 
10  since  this  is  the  only  term  affecting  the  magnitude  of  the  integrand 
and  has  a  maximum  value  equal  to  2. 

Program  ’Two’  computes  the  surface  field  for 

kr  j 

Integer  ( ■ *3+1  '  (B5) 

points  on  the  wedge  surface  or  8  points  per  wavelength  in  the  dielectric 
material. 


'Two'  reads  the  two  vector  quantities  from  the  disc  file  created  by 
'One'  and  multiplies  each  vector  by  the  appropriate  function  of  r  and 


^-jkr  sin  w 

-jkr  sin  w 
e 

n 

The  trapezoidal  rule  Is  used  to  evaluate  the  integral. 


(B6a) 

(Bbb) 

The  two  tangential 
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surface  fields  are  in  a  disc  file  as  input  to  program  'Farfld.' 

Program  'Farfld'  reads  the  two  tangential  surface  fields  from  'Two' 
and  integrates  these  vectors  in  (4.9)  using  Romberg's  method  of  integra¬ 
tion.  This  requires  integrating  the  surface  field  using  first  two  points 
per  dielectric  wavelength,  then  four  points  per  dielectric  wavelength,  and 
finally,  eight  points  per  dielectric  wavelength.  The  far  field  is  calcu¬ 
lated  for  5°  increments  from  0°  <  £  <_  180°. 
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Antenna  Parameters 


fie  slab-coordinates  and  geometry. 


wedge  antenna-coordinates  and  geometry. 


lie  wedge. 


FIG.  6. 
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FIG.  11a.  Relative  magnitude  o 
the  evanescent  region 


FIG.  12a.  Relative  magnitude  of  the  electric  field  along  the  wedge  surface  and  along  its  geometric  extension  into 
Lhe  evanescent  region  of  the  slab. 
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FIG.  13a.  Relative  magnitude  of  the  tangential  magnetic  field  along  the  wedge  surfac 


FIG.  14a.  Relative  magnitude  of  the  tangential  magnetic  field  along  the  wedge  surface  and  along  its  geometric 
extension  into  the  evanescent  region  of  the  slab. 
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of  Che  tapgential  magnetic  field  along  the  wedge  surface  and  along  its  geometric 
evanescent  region  of  the  slab. 


between  the  magnetic 


Relative  magnitude  of  the  tangential  magnetic  field  along  the 
wedge  surface  and  along  its  geometric  extension  into  the  eva¬ 
nescent  region  of  the  slab. 


difference  between  the  magnetic  field  H  and  a 
wave  along  the  wedge  surface  and  along  £ts  geo¬ 
extension  into  the  evanescent  region  of  the  si, 
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FIG.  21.  Relative  magnitude  of  the  electric  field  along  the  wedge 
surface  and  along  its  geometric  extension  into  the 
evanescent  region  of  the  slab  -  experimental  vs.  theo¬ 
retical  results. 
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de  of  the  electric  field  along  the  wedge  surface  and  along  its  geometrical  extension 
ent  region  of  the  slab  -  experimental  vs.  theoretical  results. 


FIG.  23.  Relative  magnitude  of  the  electric  field  along  the 

wedge  surface  and  along  its  geometric  extension  into 
the  evanescent  region  of  the  slab  -  experimental  vs. 
theoretical  results. 
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